Energy loss rates of hot Dirac fermions in epitaxial, exfoliated and CVD graphene 
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Energy loss rates for hot carriers in graphene have been measured using graphene produced by 
epitaxial growth on SiC, exfoliation and chemical vapour deposition (CVD). It is shown that the 
temperature dependence of the energy loss rates measured with high-field damped Shubnikov-de 
Haas oscillations, and the temperature dependence of the weak localization peak close to zero field 
correlate well, with the high-field measurements understating the energy loss rates by ~40% com- 
pared to the low-field results. The energy loss rates for all graphene samples follow a universal scaling 
of Tg at low temperatures and depend weakly on carrier density cx n~2 evidence for enhancement 
of the energy loss rate due to disorder in CVD samples. 

PACS numbers: 73.43.Qt, 72.80.Vp, 72.10.Di 



INTRODUCTION 



Despite its remarkable properties, a key problem in the 
commercialization of graphene has been that the fabri- 
cation method used for its discovery, micromechanical 
exfoliation [ij, is not amenable to large-scale commer- 
cial production. There are, however, other production 
methods which can easily be scaled up, such as epitax- 
ial growth on SiC, or chemical vapour deposition (CVD) 
on thin metal films [5]. While both these techniques are 
well established, the majority of graphene research is still 
carried out on exfoliated graphene. The principal reasons 
for this are the comparative ease of graphene fabrication 
using this method, and that to date it still produces the 
highest quality samples [3]. There is therefore a great 
need for comparisons to be made between graphene pro- 
duced by the 'research' method, and those produced in a 
more commercially amenable manner. 

Carrier energy loss rates are a particularly impor- 
tant parameter as they influence thermal dissipation 
and heat management in modern electronics as well 
as low-temperature applications such as quantum resis- 
tance metrology Q and hot-electron bolometers [11. En- 
ergy loss rates have previously been measured in ex- 
foliated graphene[B-Q, with conflicting results. There 
has also been theoretical disagreement as to how Dirac 
fermions in graphene lose energy to the lattice and how 
this varies with temperature and carrier density d 
Here we compare the temperature dependence of the car- 
rier energy loss rates in graphene produced by three dif- 
ferent fabrication methods: exfoliation, CVD and epi- 
taxial growth on SiC, for carrier densities ranging from 
1 X 10^^ cm^^ to 1.6 X 10^"^ cm^^ using two indepen- 
dent methods (weak localization (WL) and Shubnikov-de 
Haas (SdH) oscillations) and demonstrate that a single 




FIG. 1: Example R^x for the SiC/G and CVD samples. 
SiC/G, graphene grown epitaxially on SiC, before (Red, n 
= 1.63 x lO'^^cm"^, 160/.im x 35^m), after (Green, n — 4.72 x 
lO^^cm^^), and after further (Blue, n = 1 x lO^^cm^^) pho- 
tochemical gating. CVD, graphene grown by CVD onto thin 
film metal, (n = 1.62 x lO^^cm'^, 120^im x 32/im). 



consistent picture exists. 
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METHODOLOGY 



Samples 



The CVD graphene was grown on thin-film copper, 
subsequently transferred to Si/Si02[lij]- The graphene 
wafers were lithographed by e-beam into Hall-bars us- 
ing gold-only final contacting as described in our previ- 
ous work on exfoliated graphene [6]. The SiC/Graphene 
(SiC / G) was epitaxially grown on the Si-terminated face 
of SiC '12!]. Hall-bar devices were produced using e-beam 
lithography and oxygen plasma etching with large area 
titanium-gold contacting. UV exposure of top surface 



polymers was used to control the carrier density 13 1. 
The carrier mobilities as deduced from the sample re- 
sisitivity and the two dimensional carrier density were 
strongly carrier density dependent, falling from 24,000 
cm^V-^s-i at n = 1 X 10" cm^^ to 6,000 cm^V^^s-i at 
n = 1.6 X IO12 cm-2 for SiC/G and to 300 cm^V-^s-i 
at n = 1.6 X 10^'^ cm~^ for the CVD graphene. All elec- 
trical measurements were carried out using DC constant- 
current sources and multimeters. 
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FIG. 2: Kxx magnetic field traces taken in CVD graphene 
with n — 1.43 x 10^^ cm~^. All traces were taken at 500nA 
from 1.5K to lOOK. Traces are shifted such that the peak 
of the weak-localization for all traces lies at OQ. The inset 
expands the region around B = OT to show clearly the weak- 
localization peaks. 



RESULTS 
Shubnikov-deHaas oscillations 

The conventional method of measuring the tempera- 
ture dependence of the energy loss rates is by compar- 
ison of the temperature-induced damping of SdH oscil- 
lations with that induced by current heating. In both 
cases the damping is caused by thermal broadening of 
the carrier distribution. This method can be used for 
high density samples with a sufficient mobility to display 
sinusoidal oscillations. Using the well known Lifshitz- 
Kosevich formula, which has been shown to apply to both 
conventional [l^ and Dirac-like[ll] two dimensional sys- 
tems, we calculate the carrier temperature (Tg) from the 
damped amplitude with 



= Ji^^cT) . e 

p smh X 



(1) 



where B is the magnetic field and vp is the Fermi velocity, 
which is 1.1 X lO^ms"^ as measured in both epitaxial 
SiC/G[3] and exfoliated material [Tt^ . 

For the CVD material the carrier densities are nearly 
an order of magnitude higher than studied previously. 
To validate the method at these high carrier densities 
we measured the lattice temperature dependence of the 
damping, as shown in Fig. [51 using a sufficiently low cur- 
rent to avoid heating the carriers above the lattice tem- 
perature. Fig. [3] shows these amplitudes as a function of 
temperature fitted using Eq[TJ The data are well fitted 
by this equation, and the Fermi velocity, vf, is deduced 
from Eq. (3] The measured values of v^ are shown in Ta- 
blc|T]for a carrier density of 1.43 x 10^^ cm~^. The mean 
value of VF = 1-083 x 10^ ± 0.02 x 10^ ms^^ is consis- 
tent with the values reported for exfoliated and epitaxial 
graphene at lower densities which shows that the Fermi 
velocity remains essentially constant up to energies of at 
least Ef = 480meV. 

When the ambient temperature is fixed, but a current 



where r„ is the quantum lifetime and 



(2) 



and fuvc is calculated from the separation of the N and 
N-l-1 Landau levels from 



En = sgn(N) x VF\/2^hB\N'l 



(3) 



Filling 


Field 


Landau 


Velocity 


factor u 


B{T) 


level N 


Vp (ms^-'-) 


34 


17.69 


8 


1.06 X 10" 


38 


15.83 


9 


1.09 X lO** 


42 


14.33 


10 


1.10 X 10*^ 






Average vp = 


1.083 X 10" 



TABLE I: Values of vf calculated from the fitting of the tem- 
perature dependence of the Shubnikov-de Haas amplitudes, 
taken in CVD graphene with a carrier density of 1.43 x 10^^ 
cm"^. 
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P = 



nA 



(4) 



Temperature (K) 



FIG. 3: Temperature dependence of the normalized amplitude 
of Shubnikov-de Haas oscillations from a CVD sample with 
n = 1.43 X 10^'' cm~^, at the three smallest measured filling 
factors. Fits of ^.^^ ^ are shown from which the Fermi velocity 
is derived. 




FIG. 4: R^^ for CVD graphene with n = 1.43 x 10" cm'^ 
All traces taken at 1.5K with the current varied from 500nA 
up to 1.5mA. Traces are shifted such that the WL peak is at 
OQ. The inset shows the WL peak around B = OT. 



is passed through the sample, the carriers are unable to 
lose energy at a sufficient rate to reach thermal equi- 
librium with the lattice and thus heat up[l8-22[. The 
energy loss rates were determined as a function of carrier 
temperature by measuring the amplitudes of the SdH os- 
cillations as a function of current (FiglH . Comparing the 
amplitudes to the undamped value, taken at low current, 
and using Eq.([T]), the carrier temperature for each trace 
is determined. The associated energy loss rate per carrier 
(P) is calculated using 



where R^x is the longitudinal sample resistance, n is 
the carrier density, and A is the device area between the 
resistivity contacts where the power is dissipated. 



Weak Localization 

For some samples, however, the SdH method for ex- 
tracting energy loss rates could not be applied, most no- 
tably for the lower carrier density SiC/G devices, where 
the carrier densities are too low to support a series of 
sinusoidal SdH oscillations (Fig. [T|). SiC/G samples do, 
however, have prominent WL peaks. 



The magnitude of the WL peak is known [23|-|25| to 
be strongly damped by increasing temperature. En- 
ergy loss rates have previously been measured as a func- 
tion of carrier temperature from WL peaks for several 
other materials U22|. In GaAs[20| good agreement was 
found between the values for electron temperature de- 
duced from the WL peak and SdH oscillations at low 
fields (~0.5T) up to ~4K. Similarly excellent agreement 
was found for the two methods in SiGedi but could only 
be studied for electron temperatures up to I.IK. In our 
work the WL peak persists up to temperatures of almost 
lOOK allowing a much more extensive range of parame- 
ters to be studied. Previous work in graphene |26|, how- 
ever, was not able to determine energy loss rates from 
the damping of the WL peak. 

We determine the energy loss rates by comparison of 
the current and temperature dependence of the WL peak. 
The magnitudes of the peaks are measured by taking the 
difference of Rxx between OT and a fixed small field, suf- 
ficient to entirely suppress the WL behavior. The tem- 
perature dependent peak height measured at low current 
is fitted with an interpolated curve. The magnitude of 
the WL peak measured as a function of current is then 
compared with the interpolated temperature dependence 
and used to deduce the carrier temperature as a function 
of current (Fig. [S|). 

It is worth emphasizing that for the above analysis 
to be correct, the effects of raising both the carrier and 
lattice temperature (T^) need to be equivalent. For this 
to be the case the inelastic scattering time must only 
depend on the carrier temperature, and be independent 
of the lattice temperature. This has been verified by 
Ki et at by measuring the effect on of changes 
to the carrier density at temperatures up to 20K. The 
inelastic scattering, , in graphene was shown to depend 
on Coulomb interactions and Nyquist scattering, both of 
which vary with carrier temperature only. This can be 
qualitatively seen by the similar functional dependence 
upon current and temperature, shown in Fig. [S] 
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FIG. 5: An example set of data displaying the matching of the 
height of the WL peak as a function of current (blue squares) 
and temperature (red circles), from a SiC sample with n = 
4.13 X 10" cm-^ 



Energy Loss Rates 

The temperature dependence of the energy loss rates 
as measured by both techniques up to 90K are shown in 
Fig. [6] for a CVD sample which compares the results for 
the WL method and analysis using the SdH method at 
three different occupancies. This demonstrates that there 
is a power law dependence of the energy loss rates which 
tends towards T* at low temperature and that there is 
also a systematic difference between the two methods 
with the SdH technique giving energy loss rates ~40% 
lower (Table HII at any given electron temperature. This 
difference is in contrast to GaAsjio^ and SiGe^l^ where 
good agreement between the two methods was found 
when SdH oscillations below IT were analysed. We at- 
tribute this to a magnetic-field suppression of the energy 
loss rate, as the WL method measures the energy loss rate 
near OT, whereas the relatively low mobilities in most 
graphene samples require the SdH method to make mea- 
surements at much higher fields, typically lOT to 18T in 
our samples, where the Landau quantization energy is 30 
to 50 meV and significant changes in the current distri- 
bution also occur. No systematic dependence of electron 
temperature as a function of occupancy factor, could 
however, be detected at high fields, as shown in Fig. [S] 
for I'— 34, 38 and 42. Taken altogether therefore, these 
methods allow the systematic measurement of the energy 
loss rates per carrier as a function of carrier temperature 
for any graphene sample, and permit a robust compari- 
son to be made between such samples. 

Fig. [7] summarises the results from both techniques 
on all of the types of graphene, with the exception that 
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FIG. 6: Carrier energy loss rate as a function of electron tem- 
perature for a CVD sample as measured by both the WL and 
the SdH technique at occupancies of v^SA, 38 and 42. Power 
law dependences are shown for and T'^ for comparison. 



Sample Type Carrier Density SdH a WL a aSdH / 

(cm-2) (WK-4/ (WK^*/ aWL 
carrier) carrier) 

CVD 1.62 ■ 10^^ 2.4 10-^** 3.8 -10-^** 0.63 

CVD 1.43-10" 2.0 -lO-i** 3.0 -10-1* 0.67 

CVD 8.16 1012 2.2 -10-1** 3.7 10-1* 0.59 

Epitaxial 1.63 ■ lO^^ 2.1 ■IQ-i'^ 3.7 -IQ-i* 0.57 

Average = 0.62 



TABLE II: The a values P = a{Tt -Tt) for CVD and SiC/G 
samples measured by the SdH and WL methods. 
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FIG. 7: Carrier energy loss rate as a function of electron 
temperature for a representative variety of samples and from 
the two techniques used. The data follow a similar trend for 
all samples, with the energy loss rate scaling by approximately 
T*. An extrapolated T* dependence from low temperature 
calculations 9] is shown for a carrier density of 4.13 x lO^i 
cm~^ and a T'^ dependence is also shown. 
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only the SdH method could be used on the exfoliated 
graphene, which does not exhibit large WL peaksf28!]. 
Overall this demonstrates that there are no significant 
differences in the systematic behaviour between samples 
derived from the different production methods. We find 
that samples produced by all the measured production 
methods are well fitted to a low temperature, limiting 
behaviour of P = a{n){T^ — T^), where a{n) is a scaling 
constant that is weakly carrier density dependent. This 
can be approximated to the form P = a{T^) when T^, ^ 
Tg. Theoretical predictions of the energy loss rate at low 
temperatures from KubakaddiQ are shown for the exam- 
ple of a carrier density of n = 4.13 x lO^^cm"^. These 
are dominated by deformation potential coupling and are 
shown to be an excellent fit with a deformation potential 
of 19 eV. This (T^ ) power law behaviour is typical for 
low temperatures below the Bloch-Griineisen tempera- 
ture [Tbg) as observed for the resistivity in graphene |29l 
and is broadly consistent with the results of Betz et al. [8| 
who reported an approximate (Tg ) power law at high 
electron temperatures, but with a reduced coupling con- 
stant. This is in contrast to some previous theoretical ^lo| 
and experimental[3] work which suggested other powers 
of T. 

The value of Tbg {=hvskF /ks), where kp is the Fermi 
wavevector) for the samples studied is in the range 10- 
125K (1-160 X 10^^ cm^^), but as found previously in 
GaAs[l8| the data show no evidence of approaching a 
high temperature limit of linear T dependence, even 
at temperatures significantly above Tbg- At the high- 
est temperatures there is however some decrease in the 
energy loss rate, closer to a power law in the region 
of T^ or lower, as reported previously for exfoliated 
graphene [i^] at high Tg. It is likely that at higher tem- 
peratures above lOOK some additional contribution may 
occur due to optical phonons in both the graphene and 
the substrate [ifl][33|. 

Fig. [8] shows the fitted a coefficient for data below 
30K for all the measured samples, and where possible for 
multiple measurement techniques, plotted against carrier 
density. The theoretical predictions [91 of the absolute 
values are in excellent agreement with experiment for the 
epitaxial and exfoliated samples and show an a cx 
dependence on carrier density as predicted. However for 
CVD samples the energy loss rates appear to be a fac- 
tor of 2 - 3 larger. Recent theoretical work from Song 
et al.[3l| has proposed enhanced energy loss processes 
in graphene which they term "supercollisions" . These 
are caused by disorder-assisted, and two-phonon scatter- 
ing at high k values and are particularly important for 
higher energy carriers which play a disproportionately 
large role in the energy loss processes. Our CVD sam- 
ples contain significant amounts of disorder and hence the 
"supercollision" process could potentially account for the 
increased energy loss rates of the CVD samples relative to 
the theoretical predictions [il. Song et al. predict that for 
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FIG. 8: The a values as a function of carrier density for the 
CVD and epitaxial graphene data deduced from WL, and 
SdH methods. The point labelled in the key as exfoliated 
(Est. WL)is estimated from previous SdH results for an ex- 
foliated sample[^ by adjusting the value upwards by a factor 
1.6 to take account of the the measured systematic differ- 
ence between the two methods(Table [H]). Solid line on the 
graph is the theoretical prediction from Kubakaddi 9] using a 
deformation-potential coupling constant of 19 eV. 

strongly disordered samples (using kpl ~ 2) this should 
give a dependence oc T^ above Tbg, which has recently 
received experimental support from measurements of op- 
tical heating |32,] and noise thermometry [s^. Fitting the 
higher temperature data for the CVD samples in Fig. [7] 
gives P oc T'^, suggesting that this process is beginning 
to contribute to the energy loss rate. 

CONCLUSIONS 

In conclusion we have performed an extensive compar- 
ison of the transport properties of graphene produced by 
the three most common production methods. The energy 
loss rates per carrier as a function of carrier temperature 
follow a low temperature limiting T^ dependence in all 
cases, with a multiplicative factor that varies weakly with 
carrier density in good agreement with theoryj^ based 
on deformation potential scattering. The close agree- 
ment between epitaxial graphene grown on SiC and exfo- 
liated fiakes deposited on Si02 suggest that the substrate 
plays very little role in determining the energy loss rates. 
An enhanced energy loss rate is observed in CVD sam- 
ples deposited onto Si02 which may be attributable to 
"supercollisions" [3l| allowed by the presence of increased 
disorder and which may be beneficial for high power elec- 
tronic applications of graphene. We have also established 
that energy loss rates can be measured over a wide tem- 
perature range using the WL correspondence method and 
correlate well with the conventional SdH method, but the 
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values deduced at high fields (10-18T) are suppressed by 
~40% relative to the low field values. 
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